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The phase diagram of the simplest approximation to Double-Exchange systems, the bosonic 
Double-Exchange model with antiferromagnetic super-exchange coupling, is fully worked out by 
means of Monte Carlo simulations, large-A'' expansions and Variational Mean-Field calculations. 
We find a rich phase diagram, with no first-order phase transitions. The most surprising finding is 
the existence of a segment like ordered phase at low temperature for intermediate AFM coupling 
which cannot be detected m neutron-scattering experiments. This is signaled by a maximum (a 
cusp) in the specific heat. Below the phase-transition, only short-range ordering would be found 
in neutron-scattering. Researchers looking for a Quantum Critical Point in manganites should be 
wary of this possibility. Finite-Size Scaling estimates of critical exponents are presented, although 
large scaling corrections are present in the reachable lattice sizes. 

PACS numbers: 75.30.Et, 64.60.Fr, 05.10.Ln. 



I. INTRODUCTION 



There are at least two motivations for studying the 
spin-only version of the Double-Exchange (DE) models. 
On the one hand one, has its relationship with the Colos- 
sal Magnetorresistance effect (CMR)iiiSi^ On the other 
hand, in this problem some puzzles arisa^ with the Uni- 
versality Hypothesis that deserve a detailed study. Let 
us start addressing the first aspect. 

CMR has renewed the interest in Double-Exchange 
(DE) systemsi^ The typical CMR manganites are 
Lai-xAE^Mni-yOs , where ^_E'=Ca, Sr in the range 
0.2 < a; < 0.5. It is believed that the relevant degrees 
of freedom^ arc the localized S = 3/2 Mn'^+ core spins, 
and the eg holes. The Mn'^^ ions form a single cubic lat- 
tice and, besides the DE mechanism, interact through an 
antiferromagnetic (AFM) superexchange coupling. The 
relatively high spin of the Mn'^+ core suggests to treat 
them as classical spins, (f>i. Although phonons are be- 
lieved to be crucial for the CMR effect manganites 
display a very rich magnetic phase diagram which can 
be addressed neglecting lattice effectsi^ In spite of these 
simplifications, and of the introduction of powerful new 
tools)2iifliii the numerical study of the DE model in large 
lattices beyond Mean Field Approximation is out of reach 
for present day computers. Yet, finite-size effects in these 
systems are unusually large^ The need to obtain reliable 
predictions has made people to further simplify models, 
replacing eg holes by an effective interaction among the 
localized 5 = 3/2 Mn'^^ core spins. Indeed, a simple 
calculation^ shows that the kinetic energy of the elec- 
trons depends on the relative orientation of neighboring 



Mn'^"'" core spins as yl + cficfj. This substitution of a 
simpler spin-only problem in place of the very difficult 
electronic problem lies at the heart of several theoreti- 
cal analyses (see e.g. de Geniics in Ref. 0) and numer- 
ical simulations^ In spite of this, up to our knowledge 
there is only one detailed previous studjii^ of the phase- 
diagram of the bosonic DE model. That study predicted 
the existence of a disordered (PM) phase at very low 
temperatures for intermediate super-exchange coupling. 
This is very remindful of the presence of a Quantum- 
Critical Poin^ which is believed to be of importance 
for the CMR phenomenonji^ and has been predicted to 
occur in manganites by some model calculations^ The 
experimental characterization of this Quantum Critical 
Point is a wedge of paramagnetic (PM) phase, maybe 
glassyjii that at zero temperature becomes a single point 
separating two ordered phases>i^ The glassy wedge would 
be created by disorder^ and would be separated from 
the Paramagnetic state at high-temperature scale T*. 
Maybe the most surprising result of the here presented 
analysis is that this glassy wedge could not be PM or 
glassy at all, but ordered in a segment-like wajii2ii2i2£ 
(as in liquid crystals). This ordering will be referred to 
in the following as RP^ (real projective space). As we 
shall show, the RP^ phase cannot be detected in neutron- 
scattering experiments (although a short-range ordering 
will be present). Nevertheless, the phase-transition can 
be studied experimentally using the specific- heat, that 
should present a maximum (furthermore, a cusp) at the 
critical temperature. Indeed, the thermal critical ex- 
ponent is predicted^ to be = 0.78(2) which implies 
a = —0.34, and hence the cusp behavior follows. An- 
other bonus of our simplified model is that it allows us 
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to study qualitatively (see subsection IIV C|l the funny 
interplay between ferromagnetism, antiferromagnetism, 
temperature and applied magnetic field in low-doped 
Lai_,Sr,Mn03M*2i24 

Let us now address Universality. A common wisdom 
is that the critical properties of a system arc given by 
its dimensionality and the local properties (i.e. near the 
identity element) of the coset space Q /H, where G is the 
symmetry group of the Hamiltonian (the symmetry of the 
high temperature phase) and Ti. is the remaining symme- 
try group of the broken phase (low temperature). So, 
systems with locally isomorphic Q /H belong to the same 
Universality Class. This seems to be true in perturbation 
theory, where the observables are computed by doing sc- 
ries expansions around the identity element oiQ/Ti.. In 
this picture, a phase transition of a vector model, with 
0(3) global symmetry and with an 0(2) low tempera- 
ture phase symmetry, in three dimensions must belong 
to the 0(3)/0(2) scheme of symmetry breaking (classi- 
cal Heisenberg model). In addition, if 7i~0{\)~7j2 is the 
remaining symmetry, the corresponding scheme should 
be2^0(4)/0(3) which is locally isomorphic to 0(3)/0(l). 

Hence, it is interesting to check if the global properties 
of the coset space Q/T-L are relevant or not to the phase 
transition. The common wisdom has been challenged in 
the past by the so-called chiral models However, the 
situation is still hotly debated: some authors believe that 
the chiral transitions are weakly first-order while oth- 
ers claim^^ that the Chiral Universality Class exists, im- 
plying the relevance of the global properties of Q/H.. On 
the other hand, we do not have any doubt about the sec- 
ond order nature of the PM-RP^ transition. A detailed 
study of the critical exponents was recently published^ 
in Letter form. In the present work, we perform a de- 
tailed Monte Carlo, Mean Field and large-A^ study of 
the phase diagram. A thorough study is performed of 
the RP^ phase. We shall confirm that the pattern of 
symmetry breaking is 0(3)/0(2), implying a violation 
of Universality for which wc do not have any plausible 
explanation. 

The layout of the rest of this paper is as follows. In sec- 
tion ^jwc define the model, study the phase diagram at 
zero temperature, and define the order parameters and 
observables measured in the Monte Carlo simulations. 
The Mean- Field calculation is explained in section IIIII 
where we also report the results of the Large- iV analysis. 
In section Hvl we present our Monte Carlo results. We 
start determining the global phase diagram in subsec- 
tion The RP^ phase is investigated in more detail 



in subsection llV Bl while the effects of a magnetic field on 
conductivity close to a ferromagnetic-antifcrromagnetic 
transition is considered in subsection IIV CI We present 
our conclusions in section We complement the paper 
with two appendices. Appendix ^ contains the details 
about the large- calculation. In Appendix iBlthe reader 
will find the Mean-Field phase diagram as obtained from 
the fourth-order expansion of the Free-Energy. 



II. THE MODEL 



The Hamiltonian 



We define a system of spins {4>i\ living in a three di- 
mensional cubic lattice of size L (and volume V = L^) 
with periodic boundary conditions. The spins are three- 
component real unit vectors. We consider the Hamilto- 



H = - 



<4,i> 



(1) 



where the sum is extended to all pairs of nearest neigh- 
bors and we consider only J < 0. Notice that we will 
measure temperature in units of the double exchange con- 
stant. The cubic lattice is bipartite, therefore we shall 
call the lattice site i even or odd according to the parity 
of the sum of its coordinates, Xi + yi -\- Zi . 

We will consider the system at a temperature T, the 
partition function being 



Z = 



-HIT 



(2) 



where the integration measure is the standard measure 
on the unit sphere. 



B. Phase diagram at zero temperature 

As usual, the study of the phase diagram begins with 
an understanding of the ordered phases at zero tempera- 
ture. We can write in a compact way our original Hamil- 
tonian: 



<i-J> 



where 



V{v) = Jy+^l + y, 



(3) 



(4) 



and clearly y G [—1, 1]. In the limit of zero temperature, 
the only configurations which contribute to the partition 
function arc those which provide a maximum of V{y). If, 
as confirmed by the Monte Carlo simulations, the spin 
texture itself is bipartite, the value of y will be uniform 
through the lattice. Thus, a simple computation yields 
that the maxima of V{y) are at the following values of y 
(denoted by y^ax) 



1 



-1 



1 

TP 



for J> 
for J < 



1 
1 

2V2 ■ 



(5) 



It is clear that j/max — 1 corresponds to a ferromagnetic 
state and that in the J —00 limit we reach an anti- 
ferromagnetic one (j/max = ~1)- The intermediate values 
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of ymax correspond to a ferrimagnet if < ^max < 1 
and to an anti- ferrimagnet when — 1 < j/,nax < 0- The 
physical picture is as follows: the spins in the, say. even 
sublattice are all parallel along (for instance) the Z-axis. 
On the other hand the odd spins lie on a cone forming 
an angle {cosd = ymax) with the Z-axis. 
The corresponding free energy is just 



f{J) 



\/2+ J for J > --i= 



1 

2\/2 

-Tl-J for J<-^. 



(6) 



Hence we have the following phase transitions: 

1. Ferromagnetic-ferrimagnetic at J = — l/VS. It 
is easy to check that d//dJ is continuous at 
J = but d^//dJ^ is discontinuous. Hence, 
according to the standard Erhenfest classification, 
we have a second order phase transition. 

2. Ferrimagnetic-antiferrimagnetic at J = ^1/2, 
where the free-energy is C°°. At this special point 
2/max changes from positive to negative. The fact 
that ymax — implies that one can reverse ev- 
ery single spin independently of the others without 
changing the energy (more pedantically, one finds 
a dynamically generated Z2 gauge symmetry^). 

3. The limiting value ?/max — that corresponds to 
an antiferromagnet rather than an antiferrimagnet 
is reached only at J = —00. 

The transition 2 (Ferrimagnet-antiferrimagnet) needs 
further discussion. We can expand the Hamiltonian 
around the minimum y = 0, and we obtain 



1 



(7) 



Thus, close to J = —0.5 and T = one has, neglecting 
constant terms, 



O 



(8) 



<ij> 



which corresponds to an anii-ferromagnetic RP^-theory. 
The minimum energy configuration satisfies y = 0. 
Hence, we obtain that the ferrimagnet-antiferrimagnet 
transition occurs at zero temperature via a RP^ state at 
a single point. We shall see that at finite temperature the 
RP^ phase occupies a region (rather than a single point) 
of the phase diagram. 

From the previous analysis at zero temperature, one 
expects to find the following phases at finite temperature: 

• PM: the usual disordered state, where all the sym- 
metries of the model are preserved. 

• FM: a standard ferromagnetic ordering, i.e., the 
spin fluctuates around (0, 0, 1). 



• AFM: a standard anti-ferromagnetic ordering. 
Even (odd) spins fluctuate around 0° = 
(0,0,1) (<?° = (0,0,-1)). 

• FI: The ordering consists on even spins fluctuating 
around the Z axis and odd spins fluctuating around 
the cone of angle 9 < tt/2 with axis Z. 

• AFI: This ordering is similar to the previous one, 
with 9 > tt/2. 

• RP^: Here the ordering is the finite T version of the 
one found analytically inJ = — 0.5,T = 0,i.e., even 
spins fluctuating around the Z axis, with random 
sense, and odd spins fluctuating around the cone 
with random sense. 



C. Order Parameters 

In models with antiferromagnetic couplings, one might 
expect an even-odd structure of the ordered phases. 
Therefore, from the local field {4>i}, we define the stan- 
dard magnetization as the Fourier transform at momen- 
tum 0, and the staggered magnetization as the Fourier 
transform at momentum (7r,7r,7r): 



M 



Ms 



L3 



(9) 
(10) 



In a finite lattice we must take the modulus before 
taking the mean value. We will study: 



(iiM.ii) 

The associated susceptibilities are 



(11) 
(12) 

(13) 
(14) 



In order to explore RP^ type phases we introduce the 
tensor invariant under the local spin reversal. In this case 
we use as local field the matrices {ti}, constructed as 



ir'^; a,/?= 1,2,3. 



(15) 



Notice that they are traceless, thus they represent objects 
of spin 2. We can now define the associated traceless 
tensor magnetizations 

i 
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and the mean values: 

T 



^\/TrM2^ , 



The corresponding susceptibiUties are 



X 



^^(TrM^) . 



(18) 
(19) 

(20) 
(21) 



Let us close this subsection recalling the value of the 
order parameters (in the large volume limit) in each of 
the ordered phases found in the previous subsection: 

FM : > 0, = (^ fi^ > 0, = 0) , 

AFM : fiY >0, = fi'^ > 0, = 0) , 

FI : > > /.T, /ij > 0) , (22) 

AFI : > > (^ ^tT^ > 0) , 

RP^ : > A*'^ > 0, /i^ = = . 



Near a (continuous) phase transition where the corre- 
sponding correlation length, ^, diverges, the correlation 
functions in the thermodynamic limit behave for small 
p^S,^, as 



G(p) 



(31) 



Here ^ diverges as \t\~^ , t being the reduced tempera- 
ture. The anomalous dimension, t], will depend on the 
considered field. 

In a finite lattice, to estimate the correlation length 
one uses the propagator at zero momentum and at the 
minimum non-zero momentum compatible with bound- 
ary conditions. Defining F = G'(27r/L, 0, 0) and noting 
that X = G{0), one haaSS 



x/F-1 
4sin2(7r/i) 



1/2 



(32) 



III. MEAN FIELD CALCULATION 



D. Correlation Length 

For an antiferromagnetic model, not only the usual sus- 
ceptibility, also the staggered susceptibility can diverge. 
Thus, in the Brillouin zone, one needs to monitor the be- 
havior of the Green functions close to the origin as well as 
close to (tt, tt, tt). Since in critical-phenomena studies one 
usually considers only the behavior around zero momen- 
tum, we have found it convenient to define four Green 
functions in terms of four fields in momentum space 

kp) = E^"'""' ' (23) 

i 

k{p) = ^e-'P-'(-l)-»+'^-+^-,?, , (24) 

i 

f{p) = ^e-'P"-' n , (25) 

i 

fsip) - ^e-'P"''(-l)"'+^'+^'n , (26) 

i 

the Fourier transforms of the correlation functions being 



G^(p) = 


1 

u 


{kp)4*{p)), 


(27) 


GYip) = 


1 

u 


{lip)-k{p)), 


(28) 


G^ip) = 


1 

l3 


(Trt(p)tt(p)), 


(29) 


Glip) = 


1 

l3 


(Trt.(p)tt(p)). 


(30) 



Notice that GY''^{p) = G^''^{p+ (7r,7r,7r)), so that 
one could consider only non staggered correlation func- 
tions that would be studied both close to (0, 0, 0) and to 
(7r,7r,7r). 



When several phases compete, it is quite tricky to cal- 
culate the phase diagram in Mean-Field approximation 
(the T = calculation has shown that we should face this 
problem!). Since one can find different ordered phases at 
low temperatures within different Mean-Field schemes, it 
is necessary to decide which phase shall be the most sta- 
ble one. We consider that the cleanest way of performing 
such a calculation is to use the variational formulation of 
the Mean-Field approximation (see for example Ref . Isol) . 
with a variational family large enough to take into ac- 
count all the phases found in the phase diagram. In this 
way, all the phases compete on the same grounds and one 
has an objective criterion to decide which phase is to be 
found in a given region of the phase-diagram. 

One needs to compare the actual system with a sim- 
plified model where all degrees of freedom are statis- 
tically independent. The method is derived from the 
inequality^ 

F<Fo + {H- Ho)o ■ (33) 

Here, Hq is a trial Hamiltonian depending on some pa- 
rameters (the mean fields) and the average (. . . )o means 
average with the Boltzmann weight corresponding to 
Hq. The right-hand side of the inequality (|33|) is min- 
imized with respect to the free parameters in Hq, then 
used as our best estimate of the Free-energy. Thus the 
task is to generalize the standard Curie-Weiss ansatz. 
Ho ^ h ^ ■ , is the component of the local spin (j)i 
along the Z-axis), to cover all the expected orderings. 

In our case, we must use the simplest possible varia- 
tional family that permits to have different orderings in 
the even and odd sublattices: 

Ho = - Y^v^m-Y^^om- (^4) 

i even i odd 
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Notice that, as far as the calculation of the (. . .)o averages 
are concerned, all spins can be considered as statistically 
independent. Thus, the mean value of an arbitrary func- 
tion of a spin placed in (say) the odd sublattice is simply 



(odd) 




dip 



(35) 



-Vo(r)/T 



- (v/1 - (0^)2 cos ^, - (0--)2 sin^, 0^) . (36) 

We now need to parametrize the local potentials with 
the help of the mean fields, that will be our minimiz- 
ing parameters. One easily sees that keeping only the 
linear term {Vc.o{4>^) = /io,o'/'^) will not reproduce the 
ferrimagnetic or antiferrimagnetic phases, since at very 
low temperatures and non vanishing mean-fields, /ie,o, 
the spins would always be (anti) aligned with the Z axis. 
If one keeps also the quadratic term, Vc.o{4>^') = hc.o(t>^ + 
^e,o{4'^)'^, the situation improves significantly. The min- 
imum of Vc.o can now be — 1 < < 1 which implies 
that at low temperature spins would lie on the cone of 
angle 9, cos 9 = 0^,;^. Therefore, we will choose as our 
variational family 



Ho = -J2 {h,^l + U4>l?) 

i even 

- 5](Mf + Ao(0D')- 



(37) 



i odd 



As an extra-bonus, we find that the RP^-phase can be 
represented by this ansatz if the Mean Fields that mini- 
mize the r.h.s. of inequality ()33|l — at those particular T 
and J values — happen to be he = ho = 0, Ae = — Ao > 0. 
This can be explicitly checked by calculating the order 
parameters as a function of the mean-fields. Due to the 
symmetry between the even and odd sublattices, the ex- 
pressions simplify in terms of the natural linear combi- 
nations of the Mean-Fields he, ho, X^, Xo- 



h = {ho + ho)/2, 

hs = {ho-ho)/2, 

A - (Ac + Ao)/2, 

As = (Ae-Ao)/2. 



(38) 



In terms of these variables, by means of a series expansion 
in h, hg, X and As, one gets for the order parameters: 

1 

2 



.V 



/,2\ (even) 
/O 



'-Ph + ^{hX + /isAs) + 0(/l^ hi X^, Xl) 



(39) 



45 



\ (even) 
/O 



^(odd)^ 



45 



{hsX 



') (40) 

hx,) + o{hlhlxlxl), 



2 \ (even) 
/O 



((0 



3s,2\(^0dd) 



(41) 



^,P^ + l,P'(h' + h!) + ^^f3\X^ + Xl) 
+0{h\hlX\Xl), 



Ms 



((0 



= 4/3As + ^P'^hhg + 



z-i2\(odd) 




(42) 



45""'° ■ 45'' 
+0{hlhlX\Xl). 



945 



/3'AAs 



With this information in hand one can identify the differ- 
ent phases that we found at T = in terms of the non- 
vanishing Mean-Fields (of course the high-temperature 
PM phase corresponds to the vanishing of all four mean- 
fields!): 



FM : 
AFM 
FI : 
AFI : 
Rp2 : 



h. As 



hs,X > 
As > 



h = 



X 


= As 


= 0, 


X 


= As 


= 0, 


hs 


= A 


= 0, 


h 


= As 


= 0, 


hs 


= A 


= 0. 



(43) 



Let us now describe the actual calculation. As previ- 
ously said, wc introduce the function 



$(/i, hs,X,Xs)=Fo + {H-Ho)Q. 



(44) 



that, at its minimum as a hmction of ft,, /is, A and As, 
we shall identify (in Mean-Field approximation) with the 
equilibrium free-energy The partition function can be 
factorizcd to the contribution of the V/2 points of the 
even sublattice and the V/2 points of the odd sublattice: 



where 



Fo = -^(logZe+logZo). 

The average of the Mean Field Hamiltonian is 



V 
2" 



(45) 

(46) 
(47) 

(48) 



fte(r)r"nAe((0^)2)[,°™"' 
+ /io(0^)r'^+Ao((0^)2)r'^ 
As for the average of the true Hamiltonian, one finds: 

In the above expression, 0o □ is a generic spin belonging 
to the even (odd) sublattice. The problem is that, even if 
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and (/)o are statistically independent, the calculation of 
the mean- value of the square root in Eq. (|50|) cannot be 
straightforwardly factorized in an even and an odd con- 
tribution. In order to achieve this factorization, we shall 
use the series expansions introduced by de Gennesi^ One 
first use an expansion in Legendre Polynomials: 



1 



1=0 



Ai = (-1) 



i+i 



2v^ 



(50) 



(51) 



(2/- l)(2/ + 3) ' 

We can now factorize the Legendre Polynomials using 
their expression in terms of spherical harmonics 

I 



P/(0c-0o) 



47r 



21 



I E Yr{€,^.)yr(.4^i.vo). (52) 



Thus, the mean-values are factorized into even and odd 
contributions. Due to the rotational symmetry along the 
Z axis, only the to = terms in Eq. I|52|l are non vanish- 
ing. Thus we obtain. 



,(odd) 



(53) 



Fortunately, if one wants to calculate the free-energy <I> 
as a series expansion in the Mean-Fields, /i, /ig, A and 
As, at a given order only a finite number of terms in 
Eq. H5U|I contribute, due to the orthogonality properties 
of the Legendre Polynomials. This expansion allows to 
discuss the continuous phase transitions from the PM 
phase (where /i = /ig = A = As = is the absolute 
minimum of the free-energy), to ordered phases. Indeed, 
calculating $ (per unit volume) to second order one gets: 



V 



$(/i,/is,A,As) « ( 



,J0^ 
3 



2V2/32 



.P J/32 
+ {- + ^— 
^6 3 



15 
2v^ 
15 



)h^ (54) 



)hl 



+ (- 



^ 1575 
8^2/32 

1575 



^)A2 
45^ 



45 



)As^ 



This is a quadratic form in h, hs, A and As. If the above 
quadratic form is positive definite, the PM phase is a (lo- 
cal) minimum of the free energy. The other way around, 
when one of the eigenvalues of the quadratic form is neg- 
ative, the PM phase is unstable with respect to some 
ordered phase, depending on the Mean-Field that should 
grow in order to minimize the free-energy. Notice also 
that the eigenvalue corresponding to A^ is always pos- 
itive. Thus, even if there are four eigenvalues, we ob- 
tain three lines of continuous phase transitions, where 
the eigenvalues vanish: 



PM-FM line: 
PM-AFM line: 
PM-RP2 fine: 



T = 2J + 4V2/5 , 
r = -2J-4V^/5, 



(55) 



AFM 



PM 



FM 



RP^ 



FIG. 1: Phase diagram as obtained from the second order 
series expansion of the free-energy 15511 . The paramagnetic 
phase is unstable for temperatures below the full-lines (the 
instability being toward the FM, AFM or RP^ phases, as 
indicated in the plot). The dashed lines indicate the places 
where some of the eigenvalues of the quadratic form in (I55II 
vanish, but they do not correspond to phase transitions. 



Therefore the PM phase, stable at high-temperature, 
meets two transition lines of opposite slope, and an hor- 
izontal line that separates it from the RP2 phase (see 
Fig.HJ. 

For temperatures below the full lines in Fig. one 
needs to discuss the stability of a minimum of the free- 
energy different from ft, = /ig = A = As = 0. To locate 
that minimum, and to discuss its stability, one needs to 
extend the series expansion in H55|l at least to fourth- 
order in h, hs, A and Ag. This can be done (see Ap- 
pendix , but it is not particularly illuminating since 
the series expansion for $ is not fastly convergent. We 
have rather turned to a numerical method. Given a par- 
ticular value of the Mean-Fields, h, h^, A, As, we have 
calculated <& by means of a Gauss-Legendre integration 
of all the terms in Eq. (|44(l . To do this, we have divided 
the interval [— 1 , 1] in twelve subintervals and we have 
done a twelfth order Gauss-Legendre integration in each 
of them. The series of Eq. (|5()|) has been evaluated to or- 
der 50. Being able to calculate $, the minimization has 
been done using a conjugate gradient method. The re- 
sulting phase diagram is shown together with the results 
of the Monte Carlo calculation in Fig. 12 Notice that, as 
usual, the Mean-Field calculation overestimates the criti- 
cal temperatures by a large factor (about 2.3 in this case). 
However, considering this rescaling of the temperatures, 
the obtained phase-diagram is in remarkable agreement 
with the numerical one. We have confirmed that all the 
transitions are second order except the ferromagnetic- 
RP2 which is first order (nevertheless, this transition line 
is an artifact of the Mean-Field solution: in the Monte 
Carlo phase-diagram it seems to collapse to a tetracrit- 
ical point). The second order nature of the transitions 
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can be checked by computing the appropriate order pa- 
rameter at a given value of T and J, then noticing that it 
vanishes at the transition hne with Mean-Field exponents 
[M oz\T -T,\^/^ orcx |J- Je|i/2)^ 

Since Mean-Field overestimates critical temperatures, 
it is interesting to compare the previous results with the 
ones of another approximation (large N) that usually 
underestimates them. We have calculated the position of 
the PM-FM and PM-AFM phase-transition in the large 
N approximations (see Appendix fX|l : 



PM-FM Hne: T : 
PM-AFM Hne: T : 



-1.2578 J - 
-1.2578 J- 



0.5578, 
0.793. 



(56) 



The critical temperature is underestimated by roughly 
the same factor that the Mean Field approximation over- 
estimates it (see below). To extend further this calcula- 
tions would require to study non traslationally invariant 
saddle-points, which is rather complex. 



IV. MONTE CARLO SIMULATION 

The model can be investigated using a standard 
Monte Carlo method. We shall here describe some tech- 
nical points, the results being discussed in the following 
subsections. 

The single Monte Carlo (MC) step consists of a full- 
lattice Metropolis lattice sweep. Some of the simulations 
have been done at extremely low temperatures, thus the 
method of choice would have been a heat-bath algorithm, 
but its implementation in this model is rather complex. 
Fortunately, one can effectively falsify a heat-bath al- 
gorithm by means of a multi-hit Metropolis method, 
proposing per each hit as spin update a random spin 
on the unit sphere. Luckily enough, to achieve a 50% 
acceptance the number of needed hits is quite modest 
except for the lowest temperatures that represent a neg- 
ligible fraction of the total CPU time devoted to the 
problem. The pseudo-random number generator was the 
Congruential-|-Parisi-Rapuano (see e.g. Ref . . 

To extract critical exponents and critical temperatures, 
we have used the quotient methods for a pair of lat- 
tices of sizes L and 2L we choose the temperature where 
the correlation-lengths in units of the lattice-size coin- 
cide (2^L = ^2l)- Up to scaHng corrections, the match- 
ing temperature is the critical point. Let now O be 
a generic observable diverging at the critical point like 
Then, one has (up to scaling correctionsiSiS) : 

{0)2L 



{0)l 



S2L _ 



(57) 



where v is the critical exponent for the correlation length 
itself. For extracting v we have used the temperature 
derivative of the correlation-length, XQr^^ = 1 + v. To ful- 
fill the matching condition 2^^ = ^2L one often needs to 
extrapolate from the simulation temperature to a nearby 
one. This has been done using a reweighting method (see 
e.g. Ref.H. 
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FIG. 2: Phase diagram of the model. Circles correspond to 
Monte Carlo estimates. Full lines correspond to the Mean- 
Field transition temperatures divided by 2.3. The dashed ver- 
tical and horizontal lines show the direction of the Ferrenberg- 
Swendsen extrapolation in the points which are the subject 
of detailed study reported. We label to . . . ts each one of the 
more carefully studied transitions. 



A. Phase Diagram 

At a first stage, the phase diagram (see Fig.EJ has been 
explored by performing hysteresis cycles with L from 6 
to 24. The energy and relevant order parameters have 
been used to locate phase transitions. Once a rough dia- 
gram has been established, we have simulated in lattices 
L = 6, 12, and have refined the transition point calcu- 
lating the correlation lengths associated to the relevant 
order parameter and computing the point where ^/L is 
L independent. The phase diagram obtained can be seen 
in Fig. 121 We have studied in greater detail the phase 
transitions at the points to,. . . , ts, depicted in Fig. |2] 
We have simulated lattices L 6, 8, 12, 16, 24, 32, 48 and 
64, producing 20 million of MC full-lattice sweeps for 
the largest lattices in each transition. We have discarded 
5 X 10^ MC steps for thermalization. In all cases this 
has been checked to be much larger that the integrated 
autocorrelation time. In addition, at the lowest tempera- 
tures, we have compared different starting configurations 
(random, FM, etc.), concluding that the results are start- 
independent. 

Before discussing the results let us briefiy comment 
what can one expect on universality grounds. Transi- 
tion to connects the paramagnetic phase, where the full 
0(3) symmetry-group is preserved, to a FM phase where 
the symmetry group is just the 0(2) group correspond- 
ing to the global rotations around the global magneti- 
zation. Thus it is expected (and confirmed) to be in 
the Universality class of the 0(3) Non Linear a model 
(see table llll| . For all the other transitions the scheme 
of symmetry breaking is not so clear. The only obvi- 
ous symmetry-breaking (for transitions t2 and ts) is the 
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Transition 


r a 

Ij — D 


L/ — o 


r 1 o 

Lj — iz 


Lj — ID 


r OA 


r Qo 

L/ — oz 


to (T = 


0.05,^^) 










n A pi^oQi ^'9Q^ 

-u.^Oouyn^zy j 




ti [-1 = 


U.UO, fl ) 


-u.oyozoi^o J 












ti (T = 


0.05, M,^) 


-0.60083(4) 


-0.60084(3) 


-0.60078(2) 


-0.60067(1) 


-0.60052(2) 


-0.60048(2) 


t2 (J = 


-0.43, Mb") 


0.017663(12) 


0.017343(5) 


0.017163(4) 


0.017129(2) 


0.017112(2) 


0.017101(4) 


t3 (J = 


-0.8, m") 


0.07528(4) 


0.07387(2) 


0.07304(2) 


0.07283(1) 


0.07267(1) 


0.07260(1) 


t4 (T = 


0.01, m") 


-0.47199(3) 


-0.47198(2) 


-0.47196(2) 


-0.47195(1) 


-0.471919(6) 


-0.471916(3) 


t4 (T = 


0.01, Ms") 


-0.47241(3) 


-0.47219(3) 


-0.47201(2) 


-0.47196(1) 


-0.471914(6) 


-0.471912(3) 


TABLE I: Jc or Tc determined by the intersection of the correlation lengths measured in two lattices of size L and 2L. tjv(X, A) 
stands for: t]v, transition, X fixed parameter and A the order parameter associated with the correlation length considered. 


Transition 


L = 6 


L = 8 


i = 12 


L = 16 


L = 24 


L = 32 


to (T = 


0.05, m") 


0.707(4) 


0.702(7) 


0.712(12) 


0.710(10) 


0.629(95) 




ti (T = 


0.05, m") 


0.594(20) 


0.556(7) 


0.555(8) 


0.540(8) 


0.546(17) 


0.596(25) 


ti (T = 


0.05, Ms") 


0.592(5) 


0.561(6) 


0.538(5) 


0.519(7) 


0.517(13) 


0.561(22) 


t2 (J = 


-0.43, Ms") 


0.591(8) 


0.569(5) 


0.537(3) 


0.548(5) 


0.588(8) 


0.604(17) 


t3 (J = 


-0.8, m") 


0.583(10) 


0.557(4) 


0.562(3) 


0.582(6) 


0.605(7) 


0.651(20) 


t4 (T = 


0.01, m") 


0.534(4) 


0.536(10) 


0.560(10) 


0.597(15) 


0.630(17) 


0.656(24) 


t4 (T = 


0.01, Ms") 


0.545(7) 


0.564(13) 


0.581(13) 


0.611(16) 


0.629(17) 


0.650(25) 



TABLE II: Apparent v exponent obtained from the quotient method applied to {L,2L) pairs. 



symmetry between the even and odd sublattice. This is 
a Z2 symmetry, thus one might expect the transition to 
be in the Ising Universality class. The symmetries of the 
RP^ phase are intriguing and will be investigated in the 
following subsection. Let us only recall that the transi- 
tion t5 has been recently studied in great detail in Ref. 0, 
where it was found that 

Tc = 0.055895(5), (58) 
V = 0.781(18). (59) 

Perhaps not unexpectedly, the critical exponents were 
found to be compatible within errors with the one of the 
antiferromagnetic RP^ modelii^ 



Model 


V 


0(1) (Ising)*- 


0.6294(10) 


0(2)^" 


0.67155(27) 


0(3)^ 


0.710(2) 


0(4>tt 


0.749(2) 


RP^-AFMia 


0.783(11) 


Chiral(Heisenberg)— 


0.57(3) 


Chiral(XY) 


0.55(3) 


Tricriticali 


1/2 


Weak 1=* ordesiS 


1/2 


1=' order 


1/3 



TABLE III: Critical exponent v for some three dimensional 
Universality Classes. 



As for transitions tp, . . . , 14, we have located quite ac- 
curately the critical parameters (see tabled). We have 
focused in each case in the largest order parameter (see 



Eq. ^^). The PM-AFM transition should have the 
same critical behavior and we have not invested time 
in this study. We are reasonably confident in the con- 
tinuous nature of all five transitions. This stems from 
two facts. First, the energy histograms are not double- 
peaked (see an example in Fig. Yet, a much more 
refined test comes from the (L-dcpcndcnt) value of the 
effective v exponents shown in table |nj With the excep- 
tion of transition to, which as expected belongs to the 
Universality Class of the Hcisenberg model in three di- 
mensions, scaling-corrections are not even monotonous in 
their evolution with the lattice size. Although an asymp- 
totic value can not be guessed with reachable lattice sizes, 
at least one sees that, for the largest lattices, the expo- 
nent V is reasonably far from the value 1/2 to be expected 
in weak first order transitions. 



B. Detailed study of the RP^ phase 

The RP^ phase poses many questions, the first one pos- 
sibly being: is it truly a RP^ phase? In other words, (as 
far as long-distance correlations are concerned) do the 
spins behave like segments rather than arrows! As we 
will see below, the answer is positive: the spin ordering 
is truly invariant under spin reversal. Once this is es- 
tablished, one may worry about continuous symmetries. 
At the special point T = 0, J = —0.5, the spins in the 
(say) even sublattice are randomly aligned or anti-aligned 
with the (say) Z-axis. The spins in the odd sublattice 
are randomly placed on the XY-plane. This state has 
a remaining 0(2) symmetry, corresponding to rotations 
around the Z-axis. The Z2 symmetry corresponding to 
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exchanging the roles of the even and odd sublattices is 
obviously broken. Hence one may wonder if, up to the 
critical temperature which separates it from the para- 
magnetic state, the RP^ phase can be characterized as 
an 0(2)-symmetric phase with broken even-odd symme- 
try. Doubts arise from the order-from-disorder interac- 
tion that favors the alignmenti^ (in the segment sense) 
of all the spins in the planar sublattice (the odd sub- 
lattice in the above discussion). This global alignment 
would imply a second phase transition separating the low 
temperature 0(2)-symmetric phase from a phase with- 
out any remaining rotational symmetry at higher tem- 
perature. This second RP^ phase would be separated 
from the paramagnetic state where the full rotational 
symmetry, 0(3), is maintained. This complex scenario 
would explain the exotic values of the critical exponents 
corresponding to the RP^-PM transitions. Universality 
arguments^^ would then predict that the critical expo- 
nents would be the ones of the 0(4) classical Non-Linear 
a model, which are not very far from the ones actually 
foundi. On the other hand, in the simple scenario of 
a single 0(2) symmetric RP^ phase, the expected criti- 
cal exponents would be those of the classical Heisenberg 
model, in plain contradiction with our numerical results. 
Since the simple scenario seems to be the real one, we are 
facing a failure of the Universality Hypothesis for which 
we do not have any plausible explanation. 

To summarize, in the following we shall address, in this 
order, the three following questions: 

1. Is the RP^ phase truly segment-like? 

2. Is the even-odd symmetry broken up to the tem- 
perature separating the RP^ phase from the Para- 
magnetic state? 

3. Is the low temperature 0(2) symmetric RP^ phase 
preserved up to the temperature separating the 



RP^ phase from the Paramagnetic state? 



1. Tensor versus vector ordering 

We have called RP^ the phase in which the vector mag- 
netization vanishes (for any momentum in the Brillouin 
zone), and the tensor magnetization is non- vanishing, 
both at momentum (0, 0, 0) (/i"*") and at momentum 
(tt, 7r,7r) (fij). In Fig. 21 (upper and middle part) we 
show, fixing J = —0.5, that, for temperatures ranging 
from 0.001 to 0.05, there is a non vanishing thermody- 
namic limit for both quantities. For comparison, we show 
in the lower part the vector magnetization at momentum 
(tt, TT, tt) (fJ-Y), which goes to zero as We have 

also measured the correlation-length in the vector chan- 
nel. Although some short-range ordering is present, the 
correlation-length is not larger than 0.3 lattice spacing. 




FIG. 4: Lattice size dependence of fij , ^ and ^ at different 
temperatures, for J = —0.5. 

To confirm the absence of any other vectorial magne- 
tization we have measured at J = — 0.55,T = 0.5 (just 
in the middle of the RP^ phase) all Fourier components 
of the vector field 0^ for 90 statistically independent con- 
figurations for each lattice size, and plotted in Fig.^lthe 
corresponding momentum versus the maximum value of 
the Fourier component squared. In other words, we are 
searching for the maximum (over the Brillouin zone) of 
the static structure factor (divided by L"^). We have cho- 
sen as lattice sizes L = 6, 8, 12, 30, 60 to allow for different 
periodicities of the would-be vector-ordered states. If no 
vectorial ordering is present, the last quantity should go 
to zero as up to logarithmic corrections that arise 

from the fact that we are computing the maximum of a 
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FIG. 5: Scatter-plot for 90 statistically independent configu- 
rations for each lattice size {L = 6, 8, 12, 30, 60) at J = —0.55, 
T = 0.5. In the horizontal axis we plot the maximum over 
the Brillouin zone of the squared Fourier transform of the spin 
field, and in the vertical axis the corresponding associated mo- 
mentum p = \\p\\. The horizontal position of legends scales as 
Z/~^ log L in agreement with the absence-of-order prediction. 



set of 0{L^) elements. The absence of ordering is clear 
from Fig. |5l 



2. Even-Odd symmetry 

To analyze the even-odd symmetry, we measure the 
tensor correlation difference at second neighbors between 
even and odd lattices. The normalized total difference for 
a given configuration can be written as 



FIG. 6: Difference of the tensor second neighbors energies 
between sublattices, for T = 0.05 . 




3L3 



(60) 



3dd 



where the sums extend over even (odd) second neighbor 
pairs. The non-vanishing of the difference in the ther- 
modynamic limit signals even-odd symmetry breaking. 
Notice that the sublattice energy difference can be de- 
fined locally, and it plays the role of a local field. An- 
other interesting observable is the dimcnsionlcss quantity 
associated to the energy difference 



He 



(A|) 



(61) 



Fig. El shows the tensor energy difference as a function 
of temperature for several lattice sizes. A clear ther- 
modynamic limit can be observed for T < 0.05 . At 
T = 0.05 the asymptotic behavior can be elucidated by a 
direct study of the tensor energy difference histograms. A 
L = 96 lattice is necessary to clearly resolve the two peak 
structure of the histogram (sec Fig. 0) , corresponding to 



FIG. 7: Histogram of Ae for J = 0.5, T = 0.05 . 



an even-odd symmetry breaking. Notice that L = 96 
is the largest lattice used, which makes it impossible to 
study the thermodynamic limit of that quantity for T 
larger than 0.05 and less than Tc = 0.0559. We can con- 
clude that, within the computational resources employed, 
no evidence exists for a thermodynamic limit with unbro- 
ken even-odd symmetry. 

Although no thermodynamic limit can be reached 
beyond T = 0.05, more information can be obtained 
through a finite size analysis. The closer wc get to 
T = 0.05; the harder it becomes to find a two peak 
structure in the histogram. A correlation length could 
be defined in the even-odd symmetry breaking channel 
which grows as the possible critical point between the 
RP^ phase with broken even-odd symmetry and a hypo- 
thetical RP^ phase with restored even-odd symmetry is 
approached. The functional form of the growth of the 
correlation length might give indication of the existence 
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FIG. 8: Displacement of the critical temperature, defined as 
the point where ke takes a fixed value, as a function of size. 
Fits suggest that there is no even-odd symmetry restoration 
at a temperature less than the RP^-PM one. 




of such phase transition. A direct way to carry out that 
study is to define the correlation length as the lattice-size 
itself, when the histogram has a central valley at half the 
peak hight. The result shows a growth of the correlation 
length as T increases compatible with a divergence just 
at Tc, though with rather peculiar exponents. But the 
measurement of that correlation length is very noisy. A 
much more precise way to study the possible presence of 
a transition previous to Tc is to define as apparent critical 
point the T value at which ke takes a fixed value. Fig.|Sl 
shows the results. Although the possibility of an even- 
odd symmetry recovery transition previous to Tc cannot 
be discarded, results are compatible with a divergence 
just at Tc- The figures show a fit to a power law (fixing 
the critical point to the value given in Eq. H58|) 'l. It is 
worth remarking that the u exponent obtained is very 
large (2 or larger), which might point to a logarithmic 
divergence. Thus, all our results point to a single RP^ 
phase with broken even-odd symmetries at all tempera- 
tures. 



3. 0(2) Symmetry 

The chosen tool to study whether the 0(2) symme- 
try of the T = state is preserved at higher temper- 
atures, has been the eigenvalue structure of the tensor 
Mg. The latter being traceless implies that the vector 
A = (Ai, A2, A3) must lie on the x + y -\- z = Q plane. The 
whole information reduces, then, to a modulus (which is 
nothing but the observable /Xg), and an angle, which con- 
tains all the information of the eigenvalues on the sym- 
metry 0(2). As any result must be symmetric under 
eigenvalue permutations and global inversion, we can re- 
strict the angle to the interval between and 7r/6. More 



FIG. 9: Histograms of the angle of the eigenvalue vector on 
the (1,1,1) plane for two temperatures at J = —0.5. The dots 
correspond to paramagnetic configurations. 



precisely, we consider the orthonormal basis {u^, Uy\ for 
the plane given by 



- -^(-1,-1,2) . 

and define the angle 9 from the relation 



tan( 



A • Uy 

A • 



(62) 
(63) 

(64) 



with the proviso that we choose a permutation and a 
global sign such that 9 lies between and 7r/6. 

Another interesting quantity can be defined as follows. 
In the thermodynamic limit an 0(2)-symmetric phase 
corresponds to A2 = A3. We thus define 



Aa = IA2-A3I 



(65) 



which must vanish in the thermodynamic limit if the 0(2) 
symmetry is not broken. As usual the corresponding di- 
mcnsionless quantity is 



(Aa)2 



(66) 



Fig. |51 shows histograms of angles at several tempera- 
tures and lattice sizes. Dotted lines correspond to com- 
pletely disordered configurations. In case of the sys- 
tem being 0(2)-symmetric (one large eigenvalue and two 
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FIG. 10: Modulus of the difference between the two smaller 
eigenvalues of Ms as a function of T. 



identical small eigenvalues), the distribution should be 
a delta at angle tt/G. For a system with broken 0(2)- 
symmetry but unbroken even-odd-symmetry, the eigen- 
values are (a, 0, — a)-like, which would correspond to a 
Dirac delta at angle 0. Wc notice that, for small lattices, 
the distribution points to complete disorder, but as the 
size grows an inflection point turns up at T = 0.04, 0.045 
for L = 24, 48 respectively, and as L goes on growing 
a peak arises at angles ever closer to the maximum. It 
might be said, that the behavior in L is always the same, 
except for a scale change. 

Another interesting quantity is the difference between 
the two small eigenvalues (A^), which should vanish in 
the presence of 0(2) symmetry, so turning out to be an 
order parameter. Fig. 1101 shows that quantity for several 
values of the temperature and lattice size. If we look at 
an intermediate size (L = 24 for instance) the appear- 
ance is that of a transition at T ~ 0.03 to a phase with 
broken 0(2) symmetry. Yet, as L increases, the apparent 
transition moves back approaching Tc ever more. 

To check the consistency of the results with respect to 
the existence of a transition within the RP^ phase we 
can perform a Finite Size Scaling study fitting /S.\L^/'^ 
as a function of (T - To)L'^/''. Only To = Tc yields a 
reasonable fit (see Fig. lll|) . Notice that for T close to 
the definition of ceases to be meaningful, as a large 
eigenvalue exists no more since the RP^ magnetization 
fades away, and no good fit can be expected. However, 
for most T values (more precisely, for T < 0.05) the fit 
is excellent, though the t] and v values are admittedly 
rather unusual (rj = —0.5, v = 1.8). The conclusion 
should be that there is no evidence for an 0(2) breaking 
transition at any finite distance from Tc- A collapse of 
that transition over the RP^-PM transition might occur. 

A more quantitative analysis can be made studying the 
displacement of the temperature at a fixed value of k\. In 
Fig. El we plot the obtained measures together with fits 
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FIG. 11: Scaling of Ax for the analysis of a possible 0(2) 
restoring transition. Data is fairly well fitted assuming the 
transition occurs at Tc (Points next to Tc are not well fitted 
because there the largest eigenvalue becomes zero and A a 
ceases to make sense) . The fitted values are rj = —0.5, v = 1.8, 
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= 0.0559. 
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FIG. 12: Displacement of the critical temperature, defined as 
the point where nx takes a fixed value, as a function of size. 
Fits suggest that there is no 0(2) breaking transition, at a 
temperature less than the RP^-PM transition temperature. 



to several functional forms: a power law with Tc fixed to 
the value of Eq. (|58|l . a three parameter power law, and 
a Kostcrlitz-Thoulcss-likc divergence. The results point 
again to no breaking of the 0(2) symmetry inside the 
RP^ phase. 
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FIG. 13; Hysteresis along t2. The transition occurs in the 
region where f^Y changes suddenly, and the figure shows a 
movement of Tc to higher values when h = 0.01 is switched 
on (Tc(/i = 0) = 0.0171). 



C. Interplay between Ferromagnetism, 
Antiferromagnetism, Temperature and an applied 
Magnetic Field in the low-doped Lai_j;Srj;Mn03. 

In a series of papers^LSSiSi^ the interplay between 
FM, AFM, temperature and an applied magnetic field in 
the low-doped Lai-xSr^rMnOa, mainly at x close to 1/8, 
has been studied. We would like to point out some prop- 
erties of our FM-FI phase transition (point t2 in this pa- 
per) which might help to understand phenomena which, 
in those references, are related to the FM-CO (Charge 
Ordered) phase transition, not fully understood so far. 

Roughly speaking, some of the mentioned phenomena 
are: 

1. Resistance increases as T decreases below Tqo (t2 
in our model). In our simplified model, this corre- 
sponds to the fact that, when crossing the FM-FI 
transition, odd and even spins cease to be aligned, 
which makes conductivity via DEM harder. 

2. The Charge Ordered phase grows larger when an 
external magnetic field is applied. In our case, we 
have run a simulation with non zero magnetic field 
to see how the transition line moves. 

3. In the CO phase, at fixed temperature, the magne- 
tization increases with an external magnetic field, 
just as in a FM phase. 

Let us now describe the physics of the FM-FI transi- 
tion. Near the FM-FI transition, in the FM phase the 
ordering is symmetric with respect to the odd-even ex- 
change and the field fluctuates at random around the 
total magnetization vector, the fluctuations being larger 
than in the FI phase, as shown by measurements of spe- 
cific heats and susceptibilities made in both phases. More 



precisely, the magnetization increases as the temperature 
goes down from the FM to the FI phase, which can be 
explained by a diminution of fluctuations. In fact, one 
would expect that the magnetization should be smaller 
in the FI phase, with fixed odds and evens on an open 
cone around the odd direction, than in the FM phase, 
where the evens lie on a narrower cone, with a larger 
projection on the odd direction. Yet, the large fluctua- 
tions in the FM phase destroy the even contribution to 
the magnetization. The FI vacuum consists then on the 
odd, practically frozen, sublattice, and the even sublat- 
tice, with spins on an open, but less fluctuating cone. 

Let us now look in the FM phase close to the FI tran- 
sition, and switch on a weak magnetic field in the Z di- 
rection. This will have the general effect of collimating 
the spins. In more detail, odd spins shall freeze closer to 
the Z direction, which shall cause the even sublattice to 
freeze on the cone, with smaller fluctuations. Paradox- 
ically, the collimating effect of the magnetic field in the 
Z direction is to stabilize the cone, effectively opening it, 
giving rise to a more Fl-like ordering, i.e., the FI phase 
shall invade the FM phase, and the critical temperature 
shall rise. This phenomenon (see point 2) is accompanied 
by an increase in the magnetization at fixed temperature 
in the FI phase (see point 3). 

In order to check the correctness of the description, we 
have simulated in the neighborhood of the transition with 
h = 0.01, which does not alter the system properties, and 
have run a hysteresis cycle at J = —0.43 in L = 12, 24 
between T = 0.01 and T = 0.025 (i.e. along ta). A good 
observable for the transition is ji^ ■ The results at the two 
L values show that the finite size effects are negligible in 
front of the change in Tc with h. 

Fig. ^] shows the result, which confirms that the in- 
clusion of a magnetic field rises Tc, causing the invasion 
of the FI ordering into regions which at /i = were FM. 



V. CONCLUSIONS 

We have studied a simple model for double exchange 
interactions which retains a good number of interesting 
properties. It exhibits a complex phase diagram with 
ferromagnetic, ferrimagnetic phases, with their staggered 
counterparts, and a segment-ordered phase. 

We obtain quantitatively all phases with approximate 
calculations (Mean Field and 1/A^ expansions), which 
can be contrasted with exact Monte Carlo calculations. 
With Monte Carlo simulations we obtain, in addition to 
the precise positions of the transitions, information about 
their order. Our conclusion is that all transitions seem 
to be second order, although an accurate determination 
of the critical exponents is difficult and it is beyond the 
scope of this paper. 

We have studied in detail the exotic RP^ phase 
(segment-ordered), concluding that it is a single phase 
up to the resolution allowed by the lattice sizes used in 
the simulation. The presence of a RP^ phase up to T = 
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is interesting from the experimental point of view, since 
it can be confused with a PM or glassy phase and conse- 
quently with a Quantum Critical Point. We have shown 
that the structure factor {V, in fig-EJ remains bounded 
in the full Brillouin-zone . Therefore the RP^ phase can- 
not be detected in neutron scattering experiments as a 
long-range ordering, although the phase transition will 
show up as a maximum (more precisely, a cusp) of the 
specific heat. A short-range ordering would of course 
always be present. Since the critical exponent a is neg- 
ative, the Harris criterion^ implies that our results are 
robust against disorder effects. The RP^ phase is char- 
acterized by a breakdown of the even-odd symmetry and 
a remaining 0(2) symmetry, despite the fact that the 
measured critical exponents for the RP^ transition point 
to a total breakdown of the initial 0(3) symmetry. We 
have also discussed the effects of a magnetic field on the 
ferromagnetic-ferrimagnetic transition, and we have dis- 
cussed its interplay with electrical conductivity. 
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APPENDIX A: LARGE N APPROXIMATION 



with A, that action, as 

/5 



(^Ajj + Xij (fi ■ (fj - Xij pij + 2W{ptj) 



(A4) 



As we are interested (in this part of the calculation) 
in paramagnetic or/and ferromagnetic phases, we sepa- 
rate the spin in two pieces: the first one parallel to the 
symmetry breaking direction, (one degree of freedom), 
and the orthogonal part {N — 1 degrees of freedom) , (j)-^ . 
At this point, the spins have no definite modulus, and we 
can perform the functional integration over the orthogo- 
nal part of spins (a Gaussian integral) 



d[0^] 



cx exp 



N- 1 



Trlogi? , (A5) 



where R'^^ is the propagator (a, 6 = 1, A^ — 1 and i lives 
in the three-dimensional lattice) and it is given by 



rtab 



(A6) 



The sum runs back and forth along the 3 lattice axes 
and is the neighbor of site i in the direction defined 
by V. The trace, Tr, is over the space and spin compo- 
nents. The quantity \{<j)'^ ■ R4>'^) is the contribution to 
A involving the orthogonal part of the spins (which is a 
quadratic form with matrix R). 
In momentum space, R reads: 



V 



(q-q')-ri 



(A7) 



We write the model as 

'H = -N W{\ 



(Al) 



<i-J> 



The Boltzmann weight is exp(— 7i) and 
W{x) ^ Jx + ^fx . 



(A2) 



Using the expression of the Dirac deltas (one to fix the 
spin modulus, (pf = 1, and another to write that x = 

■ 4>j) in terms of functional integrals, we can write the 
partition function of the model in the following way42i^ 



Z (X / d[p, A,/i,0] e 



NA 



(A3) 



where v — — r^j. In the large N technique wc must 
maximize A. In order to keep the computation at its 
simplest level, we make an ansatz over the fields Ay , /i^, 
Pij, and (/)ll: we are assuming that we will describe under 
this ansatz translational invariant phases, like paramag- 
netic and ferromagnetic ones. So, we will consider that 
all these fields are independent of x and p and we will 
write them as A, p, p and a. Therefore, A is 



P) 
Aq log 



\cj^ + 2W{p)\+^p{l-a') 



- A ^ cos qi, 



where d is the dimension of space and 



dg = 



= 1 



(A8) 



(A9) 
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Hence, this computation is valid in paramag- 
netic/ferromagnetic phases where we have translational 
invariance. As usual we write 



p2 = 4^sin2(p,/2) . 



(AlO) 



The continuum limit of fp' is p^, and so, we can define a 
mass toq: 



2 2/i 

ttiq = — 2d . 

A 



(All) 



and A can be written as 



^ = p) + Xa^ + 2W{p)]+^ji{l-a^) 

- i log A - i I dg log [ml + f] . (A12) 

The saddle point equations are: 

!3d{l-p) + \ [{ml + 2d)I{ml) -l]+d[3(j'^ = 0, (A13) 
A 



/3(l-a2) = -/K) 



2W'{p) = A , 



(A14) 



(A15) 



longer the meaning of a mass (hence, in this case, mo = 
is not a signature of criticality). In this case the complete 
solution is: 



1 



and 



T = 



2/(0)d ' 



J 



(A20) 



(A21) 



This last equation tells us what is the magnetization 1 — 
in a given point {T,J) . In the interval J > —1/2 
we obtain the solution a ~ I. In addition in the interval 
J e (-1/2,-1/(2^2)) a second solution with a < 1 
appears. This is the signature of the ferrimagnetic phase. 
Hence, we have recovered part of the previous T = 
results. 

As mentioned, above the previous calculation is valid 
only in paramagnetic/ferromagnetic phases. In order to 
manage the paramagnetic/antiferromagnetic phases we 
use the following trick: we change the sign of the odd 
spins and we leave unchanged the even spins, so, the 
Hamiltonian reads: 



(A22) 



<jj> 



and following the technique outlined above, we obtain 
the equations of the saddle point: 



a{d\ — /i) = , 



(A16) 



A = -2W'{p) 



(A23) 



where 



I{ml) 



Aq 



(A17) 



(A24) 



One solution is a = 0, the paramagnetic phase. We 
can find a second order phase transition by fixing the 
mass mo to zero: 



dp^d{l~a^)-^{{2d + ml)I{Q)-l) , (A25) 
A 



where 



T/(0) = J 



Po 



1 



2dI{Q) 



(A18) 



(A19) 



In three dimensions /(O) ~ 0.2527. So we have found the 
critical line between the paramagnetic and ferromagnetic 
phases. This solution is only valid in J > —1/2. It is easy 
to check that the J = —1/2 vertical line corresponds to 
infinite mass. So, with these formulas we cannot reach 
the region to the left of J = —1/2. Below we will see how 
to solve this drawback. 

We can try to connect this calculation with the T = 
results. The solution cr 7^ implies that d\ = /i and so 
ttiq = 0. Notice that in a magnetized phase, mo has no 



a{d\ - /i) = 



(A26) 



Again mp = 2/i/A — 2d. 

In the paramagnetic phase ct = is the solution and 
the equation of the critical line is (obtained by fixing 



0): 



-T/(0) = J 



where 



Po = 



2dI{Q) 



(A27) 



(A28) 



In addition < cr < 1 is also a solution and so rfA = /i 
and this implies, as in the PM/FM computation, that 
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the magnetization 1 — cr^ in a given point (T, J). Again, 
in this part of the calculation we cannot reach the region 
J > -1/2. The line J = -1/2 has again = 0. 

Finally, we report the transition lines in terms of the 
temperatm'e measured in the Monte-Carlo simulation. 
Taking into account that Tmc = T/N, where T is the 
temperature of the large N calculation and fixing N to 
3, we obtain the FM-PM line: 

Tmc = 1-2578 J + 0.5578 , (A31) 

and the AFM-PM hue 

Tmc = -1-2578 J- 0.793 . (A32) 



FIG. 14: Mean-Field phase diagram as obtained from the 
numerical minimization of the free-energy (dashed lines) and 
from the minimization of the free-energy calculated to fourth 
order (full lines) . The dotted lines are artifacts of the fourth- 
order approximation. 



TOq = The phase being a (staggered) magnetized one, 
this does not imply criticality. The solution is then: 

. 1-(t2 



APPENDIX B: 



2/(0)d 



and 



T = -- 



m 



J 



2^^o^ 



(A29) 



(A30) 



As in the PM-FM case, this last equation tells us what is 



MEAN FIELD FOURTH ORDER 
ANALYSIS 



We have extended our Mean Field power expansion 
analysis to fourth order, so that we can find transitions 
where the paramagnetic phase is not involved. The an- 
alytical minimization with respect to all fields is a very 
hard task. But we can face the problem restricting the 
parameter region, using the essential fields that can de- 
scribe the transition. First of all, let us explore the transi- 
tions inside the ferromagnetic region found in the second 
order analysis. 



$/.A.(As) = $(/i""",0,0,A,), 



(Bl) 



where /i™™ is the value of h where <^h{h) = 0,0,0) 
reaches the minimum. We can expand 



^hx. (As) = 0, 0, 0) + a,a. (r, J)\l + 6,^. (T, J)\t + 0{\t) . (B2) 

Then, if b^x^ (T, J) is positive there is a stable minimum with non zero As when ah\^ (T, J) is negative. Therefore, we 
find a transition line when ah\^(T, J) = 0. In this case bh\^{T^ J) > if T < 0.31 and the transition line between the 
ferromagnetic and a ferrimagnetic phase, where M and Mg are non zero, is 

^ _ 4 (20 + 83^/2J+140J^) 

Jfm-fi — — , ^ [Ho) 

386V2 + 875 J + V369392 + 971810\/2 J + 1265425 J2 

We can do a similar analysis inside the RP^ phase. In this case wc study 

^xAh)^Hh,0,OATn, (B4) 

and 

$A.^.(/is) = *(0,/is,0,Af"), (B5) 

obtaining the following transition lines 

r 32 (327 + 406^/2J) 



35 ^480^2 + 539 J + \/-575136 - 768768^2 J + 290521J2) 
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RP2-AFI 



32 (283 + 406^2 J) 



35 (^-128V2 + 539J - ^929312 + 1148224\^J + 290521J2j 
Finally, inside the antiferromagnetic phase we find a transition to an antifcrrimagnetic ordering minimizing 



(B7) 



The transition line is 



FM-FI 



$/..A.(A.)=<i>(0,C'",0,A,) 



4 (404 + 795^/2 J + 700 



5 ( 296\/2 + 875 J + ^463688 + 1085630^/2 J + 1265425 J2) 



(B8) 



(B9) 



The fourth-order phase-diagram is depicted in Fig. ^1 together with the numerical calculation of Section. IIIII 
Letting aside the FM-RP^ line (which is first-order in Mean-Field approximation), the results of the fourth-order 
approximation are qualitative satisfying. 
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